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ABSTRACT 


A radiation shield optimization procedure based on the computer code, OPEX-n, is 
described. The OPEX-n code, based on an earlier, steepest -descent -method code 
OPEX, has been recoded to improve coding, simplify data input, use spherical geom- 
etry, and alter the dose-thickness relation when a layer has been removed. A complete 
description of how to obtain the necessary input data for OPEX-n from other transport 
calculations is given. Data input instructions, FORTRAN IV code listing, and a sample 
problem optimizing a seven-layer shield of tungsten and lithium hydride for a space 
power reactor are given. 



OP EX- II, A RADIATION SHIELD OPTIMIZATION CODE 
by Gerald P. Lahti 
Lewis Research Center 

SUMMARY 

A radiation shield optimization procedure based on the computer code, OPEX-II, is 
described. The OPEX-II code, based on an earlier, steepest -descent -method code 
OPEN, has been recoded to improve coding, simplify data input, use spherical geometry, 
and alter the dose -thickness relation when a layer has been removed. A complete de- 
scription of how to obtain the necessary input data for OPEX-II from other transport 
calculations is given. Data input instructions, FORTRAN IV code listing, and a sample 
problem optimizing a seven layer shield of tungsten and lithium hydride for a space 
power reactor are given. 


INTRODUCTION 

The radiation shield designer is faced with the task of selecting shield materials and 
material arrangements which, will not only provide adequate protection against radiation 
but will also minimize shield weight, cost, or space. If components of the total dose are 
independent of one another (i.e. , primary gammas, and fission neutrons), and the geom- 
etry is simple, analytic expressions for the radiation dose may be defined and a closed 
form solution for minimum weight may be obtained. Reference 1 reviews such cases. 

I 

When radiation groups are not independent, which is generally the case for second- 
ary gammas generated by neutron absorptions and inelastic scatters throughout the 
shield, closed form solutions are no longer attainable so numerical iterative methods 
must be employed. Two of these numerical, iterative optimization techniques are 
(1) the method of Lagrange multipliers (as applied to the shield weight optimization 
problem in ref. 2) and (2) the method of steepest descent (as applied to the shield weight 
optimization problem in ref. 3). 

In both methods, an empirical analytical expression (hereinafter called the dose- 
thickness relation) is assumed which relates the radiation dose at some reference de- 



tector point to all thicknesses of material present. The parameters in this empirical 
expression are obtained by fitting them to some accurate detailed radiation transport cal- 
culations of dose for a given base configuration and perturbations of that configuration. 
The geometry and thicknesses of material determine weight and derivative of weight with 
respect to thickness. With first derivative of weight and dose with respect to thickness 
as determined from the dose -thickness relation, the optimization procedure alters the 
base configuration to obtain a set of thicknesses corresponding to a minimum weight con- 
figuration (or at least a local minimum) for some dose constraint. Optimization codes do 
not select materials but merely alter initial configurations. Optimization codes may 
eliminate layers but they never add layers. The optimum weight estimate is only as good 
as the parameters which describe changes of dose with thickness. A final, detailed 
proof calculation is necessary to confirm the predictions of the optimization code. 

The steepest -descent method of reference 3 was incorporated in a rudimentary com- 
puter program called OPEX (ref. 4) but was limited to slab geometry. This report 
describes a revision of OPEX, called OPEX-n, and how it is applied to a'radiation 
shield optimization. The basic steepest -descent method of reference 3 has been main- 
tained, but the code OPEX has been completely rewritten to improve coding, simplify 
data input, use spherical geometry, expand the output, and alter the dose-thickness 
relation when a layer has been removed by the optimization code. In this report, a com- 
plete description of how to obtain the necessary input data for OPEX -II from other trans- 
port calculations is given. Data input instructions, FORTRAN IV code listing, and a 
sample problem optimizing a seven -layer shield of tungsten and lithium hydride for a 
space power reactor are given. 


DOSE-THICKNESS RELATION 

The total radiation dose rate D at some reference point in space is defined, for 
purposes of the optimization procedure, to be 


where 


D = 


IMAX 


D i 


D i 


IMAX 


total dose rate 

11 . 

i component of total dose rate (e. g. , dose due to capture gammas from first 
shield layer or dose due to inelastic gammas from last shield layer, etc. ) 

number of components of dose 
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Each dose component is further assumed to be of the form 

/NREG 

D i = c i ex p ^ S ^ij 

where 

Cj fitted parameter 

tj thickness of region 

l u... ’’attenuation coefficient” which describes effect of change in thickness t. on D. 

J 

NREG number of regions present 

■Hi 

The i in dose component need not be associated with the j region. There may be 
more than one contribution to the total dose from a given layer. Other layers may con- 
tribute negligible secondary gamma dose and their dose contribution is omitted. Gen- 
erally because of differences in formation rate and gamma attenuation, a given high-Z 
shield region will have one dose component due to capture gammas and a second due to 
inelastic gammas. The core will have primary neutrons and primary gammas. A hy- 
drogeneous layer may have negligible dose. 

For example, consider the spherical reactor -shield described in the following 
sketch and table with a dose constraint at point P: 
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Region j 

Description 

Radiation 
dose com- 
ponent i 

Description 

1 

Reactor core 

1 

Core neutron 



2 

Core gamma 

2 

Reflector 

3 

Capture gamma 



4 

Inelastic gamma 

3 

Lithium hydride 

— 

Negligible dose 

4 

Tungsten shield 

5 

Capture gamma 



6 

Inelastic gamma 

5 

Lithium hydride 

— 

Negligible dose 


In this example, p ^ 45 Mg 4 , and P 4 4 represent the attenuation of gamma rays by 
the tungsten layer for the various gamma sources in the core and layers between the core 
and the tungsten. The coefficient p^ 4 represents the attenuation by the tungsten of 
neutron dose due to neutrons bom in the core. In contrast, p^ ^ represents the attenua- 
tion by lithium hydride (layer 3) of neutrons which give rise to secondary gamma sources 
in the tungsten (layer 4). 

For coefficients such as M 5 4 , describing the effect of a thickness on its own source 
strength, equation ( 2 ) becomes inadequate particularly when the thickness of the layer 
goes to zero. OPEX-13 does, however, set Cj = 0 if a particular layer is eliminated by 
the code. At this point, however, the attenuation coefficients P- should be recalculated 
for the new configuration. Regions of origin of the various dose components are required 
data input to the code to facilitate this operation. 


FITTING PARAMETERS TO DOSE-THICKNESS RELATION 

The coefficients of equation (2) are obtained as precisely as possible by performing 
a series of transport calculations. A starting, base configuration is selected, prefer- 
ably as close to an optimum configuration as possible based on past experience. With 
the transport calculation, such as a discrete ordinates calculation, the individual dose 
components D. are evaluated. Each layer of the base configuration to be altered by the 
optimization code is then systematically increased by a nominal amount (say, 1 cm) and 
the dose components are reevaluated with the transport code. From equation (2), then, 
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D i^l» *2’ ' ' * V ' ' ' ’ *NREg) “ C i exp l 


D i^l» t 2 > • ' ‘ *j + At j’ • ' •’ ^NREG^ “ C i exp ( ~ ^ j p ijM ex P( _ ^ij At j) 


NREG 

j=i 

NREG 
i=i 






Solving the above pair of equations for results in 


M ti = — In 
3 Atj 


Vfi V t 2 . 


•’ v • • 


^REG^ 


D i^l’ t 2’ 


•’ V At r • 


*NREG^ 


( 3 ) 


The coefficients are determined for all i and j in this manner. 

The basic configuration data, that is, the set of thicknesses t^, dose components Dj, 
and "attenuation coefficients" constitute the required input data. For the set of 
base configuration, then, the coefficient is calculated by the code OPEX-II from 
equation (2); that is, 


/NREG \ 

C t = D i exp(2p Mijtjj 

WEIGHT-MINIMIZATION PROCEDURE 

The procedure for obtaining the minimum weight configuration by the steepest - 
descent method is presented in this section. The equations are basically from refer- 
ences 3 and 4. The narrative is expanded to illustrate the flow of computation and to 
comment on code output. 

The mathematical problem to be solved is that of minimizing the weight w, a func- 
tion of thicknesses tj, while constraining the total dose D to some particular value; 
that is, 


Minimize w(t^, t^ j • • • > ^NREG^ 


with constraints (a) D = ^ ^ C. exp 



= constant 
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(c) = constant for any desired values of & 


Constraint (b) ensures a physical solution. The optional constraint (c) is useful if it 
is desired that some thicknesses be kept from changing during the computation (e. g. , the 
reactor core and reflector thicknesses). Constraint (c) is necessary for the spherical 
geometry programmed presently to prevent the trivial case of reducing reactor core size 
to zero. 

An n-dimensional Euclidian vector space with Cartesian coordinates tj, t 2 , . . . 
t N REG defined. The following vectors are defined on this space: 


t = (tj 

*2’ * 

• • t NREG^ 


3w. 

3w 

w 

~*2 ' 

» • 5 

0t NREG, 

/ 3D 

3D 

3D 


3t 2 

• • y 

0t NREG 


The notation t = (t^, t 2 , . . . t NREG ) means tjXj +t 2 X 2 + . . . where Xj are 



unit vectors in the i direction. Vectors g and a represent the gradient of weight 
and dose. The components of g are evaluated from analytic expressions of weight as a 
function of thickness, and depend on geometry. The components of a are evaluated from 
the partial derivatives of (1) and (2), namely, 


IMAX a-n IMAX 
3D — v 0D i > 3 


0t k i=l 0t k i=l 0t k 


/NREG 


: i ex P " 


\ j=l 


li. .t. 
i] 3 


IMAX 

i=l 


NREG 


(-M ik ) c i expl- >> ] M„t. 


j=l 


i] 3 


IMAX 


^iiPl 


i=l 


The unit vector u (see ref. 3 for derivation) 




points in the direction of greatest weight decrease (steepest descent) along a hyperplane 
tangent to the hypersurface described by the equation 

D(X) = D(tp t 2 , . . . t NREG ) = Constant 

a 

Components of u, namely u., represent increments of thickness to be added to 
each tj to approach the minimum weight criterion. 

The optimization code proceeds as follows: 

(1) A fraction f of each component u- of u is added to each thickness t-. The 

J J 

fraction f is an input parameter. (A value of f = 1.0 has given satisfactory results.) 

(2) The new set of thicknesses generally does not return the correct dose constraint 
so a first-order correction is applied to each t^ to return the dose constraint. This 
correction is 

^new = ^old + [d (constraint) - D(calc)] x — 

(a • a) 

Steps 1 and 2 are repeated until the relative change in weight from one iteration to 
the next is less than some prescribed value. 

The code output includes final thicknesses and individual dose components as cal- 
culated from the dose-thickness relation for each iteration. It is incumbent on the user 
to make a final detailed proof calculation to verify the results of the prediction of the 
optimization code. Experiences have indicated that if input coefficients are determined 
accurately, final predictions of the optimization code are quite good provided the con- 
figuration is not radically changed If the configuration is changed severely, a recal- 
culation of coefficients is in order. 

SOME OTHER USES OF OPTIMIZATION CODE 

Once an optimized base case is obtained, effects of nominal changes in reactor size, 
power level, and dose constraint on shield weight may be estimated using the optimiza- 
tion code. This is done by altering reactor radius, by scaling dose components pro- 
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portionately, or by specifying a different dose constraint, respectively, and allowing 
OPEX-II to seek a new minimum weight configuration. 

Cost optimization (minimization) may also be performed by specifying cost per unit 
volume rather than density (weight per unit volume) for each region. 

Input information to the code requires knowledge of each component of the dose. Out- 
put information also includes contribution of each component to the total dose. Thus one 
learns which regions are important and most sensitive to the calculation. 


OPEX-II CODE 

In this section the code details are presented, including 

(1) Flow chart for data input 

(2) FORTRAN IV listing 

(3) Sample problem and sample problem output 

The geometry presently coded is spherical and is illustrated in sketch (b) 




The data input consists of a set of thicknesses for each region. Radii are calculated 
internally. A thickness must be specified for each region. 
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Flow Chart for Data Input 




1 


Control card 


Card 

column 

Variable 

FORMAT (315, 5F 10. 4) 

1-5 

NREG 

Number of regions (thicknesses) in this problem 
(—25) 

6-10 

MAX 

Maximum number of iterations allowed 

11-15 

IMAX 

Number of items in dose equation (<25) 

16-25 

DES 

Desired dose rate (used only if > 0.0) 

26-35 

EPS 

Convergence criterion for weight 

36-45 

EPSD 

Convergence criterion for initial dose (EPS, EPSD 
typically 0.001) 

46-55 

C0N 

Fractional step size on u (0. 5 < CON < 1.0) 

56-65 

CA 

Fractional step size for initial dose (0. 5 < CA 
< 1.0) 



T(J) Thickness of j th region; 
NREG values required. 


RH0(J) Density of j th region; 
NREG values required. 


FORMAT (7E10.4) 


FORMAT (7E10.4) 


4 


NB(J) Thickness constraint flag 
NB(J) = 0 constrain j** 1 region 
to constant thickness 
NB(J) = 1 allow j th region to 
change; NREG values required. 


FORMAT (2511) 
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i 


FORMAT (7E10.4) 


For each region j read a new card 
(or set of) cards with IMAX values 
of /LUj- NREG cards (or sets of 
cards) required. 


6 


D(i) 

i** 1 dose component corresponding to 
initial geometry read in above; 
IMAX values required. 



NU(I) 

NU(I) is the region number physically 
associated with the i dose term. 
This cross referencing is necessary 
to zero the correct C(I) in the event 
the j th region is diminished to zero 
thickness by the optimization process 
A region may be identified more than 
once. IMAX values required. 


FORMAT (7E10.4) 


FORMAT (1415) 
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FORTRAN IV Listing 


Presented in this section is the IBM-7094-11 FORTRAN IV listing for the OPEX-II 
code. SUBROUTINE WEIGHT is coded for spherical geometry. Indication of how to 
alter this subroutine for the geometry of slabs bounded by a conical surface is given in 
that subroutine. The code can handle up to 25 regions and dose components as presently 
dimensioned. Running times are typically less than 0. 1 minute on the IBM 7094-11 for 
12 -region, 16 -dose component problems. 


Sample Problem 

A sample problem consisting of a reactor, molybdenum reflector, and a shield con- 
sisting of seven alternating layers of lithium hydride and tungsten is illustrated in fig- 
ure 1. Region descriptions, densities, and initial thicknesses (guessed) are listed in 
table I. The configuration is to satisfy a dose constraint of 2 mrem per hour at a point 
20 meters distant. A series of discrete ordinates calculations for both neutrons and 
gammas was made to calculate doses from each source; perturbations were made to 
determine the attenuation coefficients p... The results of these calculations are listed 
in tables n and IE. Because the core radius, plenum, pressure vessel, and reflector 
thicknesses are to be constrained in this calculation, a value of P- = 0.0 is arbitrarily 
assigned to these regions. 

The complete computer output for this OPEX-II calculation to seek a minimum 
weight for this seven-layer configuration is given in table IV. The output consists, first, 
of a listing of all input information, followed by the value of the dose constraint (DES), 
the value of the calculated dose for the initial configuration (DOS), and the weight (WT) 
in grams of the initial configuration. If DES and DOS do not agree to within the param- 
eter EPSD, a new set of thicknesses is calculated and printed; this new set of thicknesses 
satisfies the dose constraint. The listing of DES, DOS, and WT is followed by the values 
of the dose components D(I) and the thickness of each region T(I). 

The results of each OPEX-II iteration for the present problem are shown in figure 2. 
Shown, for each iteration, is the size and relative position of each of the layers as 
adjusted in the course of the calculation. The final thicknesses and values of the dose 
components are listed in tables I and n for comparison with initial values. The initial, 
guessed configuration weighed 3. 594x10^ grams (79 000 lb); the final, 2. 999x10^ grams 
(66 000 lb). 


( 


11 



TABLE I. - REGION DESCRIPTION 


Region j 

Description 

Density, 

g/cm 3 

Thickness, 

cm 

(initial guess) 

Final prediction 
of shield 
thicknesses, 
cm 

1 

Reactor core 

9.957 

26.0 radius 


2 

Plenum 

8.647 

2.50 


3 

Pressure vessel 

16. 763 

.60 


4 

Molybdenum reflector 

9.234 

11.00 


5 

Lithium hydride 

.75 

(17.90) 

20. 52 

6 

Tungsten 

19.3 . 

(7.00) 

9. 71 

7 

Lithium hydride 

.75 

(14.00) 

12. 32 

8 

Tungsten 

19.3 

(5.00) 

2.82 

9 

Lithium hydride 

. 75 

(10.00) 

10.32 

10 

Tungsten 

19. 3 

(3.50) 

2. 33 

11 

Lithium hydride 

.75 

(59.50) 

39.29 


TABLE H. - DOSE COMPONENTS 


Region j 

Dose component, 
D i 

Value of dose component, 
Df, 

mrem/hr 

With initial 
shield 
thicknesses 

With final 
shield 
thicknesses 

1 

Neutron 

0.02430 

0.2527 

2 

Core gamma 

.00303 

.0079 

3 

Plenum, pressure vessel capture gamma 

.00196 

.0049 

4 

Plenum, pressure vessel inelastic gamma 

. 00220 

.0055 

5 

Reflector capture gamma 

.204 

.4782 

6 

Reflector inelastic gamma 

. 00504 

.0129 

7 

Region 6 tungsten capture gamma 

.0921 

.4007 

8 

Region 6 tungsten inelastic gamma 

.00974 

.0949 

9 

Region 8 tungsten capture gamma 

.0988 

.2696 

10 

Region 8 tungsten inelastic gamma 

.0278 

.0759 

11 

Region 10 tungsten capture gamma 

.201 

.2650 

12 

Region 10 tungsten inelastic gamma 

.0947 

. 1320 

Total 

0.7647 

2.000 





TABLE m. - COEFFICIENTS 
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TABLE XV. - OPEX-H LISTING AND OUTPUT FOR SAMPLE PROBLEM 


JIBFTC npt DECK 


C MAIN PROGRAM ADOPTED FROM OPEX I AI > DIST. BY RSIC/ORNL 

C OPFX BASIC EQUATIONS COMPLETELY Kh' PROGRAMMED 8 Y 


c 


G.P. 

L AH7 I 

0*NASA-LtWl S RESEARCH 

CENTER 

c 

r 




CLEVELAND , 

, OHIO 


t 

COMMON 

NR EG. 

MAX. 

IMAX. DES. 

CUN, CA, 

EPS, EPSD 


1 COS. AA. 

Wf . 

GG. 

AG. IT. 

NBI25) , 

• N I J 1 25 ) , 


2 C< 25), 

AI 25 ) 

, 

C 125) . 

E MU I 25 ,25) 

, 


3 T I 25) . 

RHOt 25) . 

G 1 2 5 ) . 

U(25> 



777 7 CALL INPUT 







CALL WEIGHT 
CALI. DOSE 

I FI D6S .GT. 0.0 ) CALL 1NIT 
DS = DOS 

DO 120 IT =1. MAX 
W = WT 

SO = SOKTIGG -AG*AG/AA ) 

A GA= AG/AA 

DO 60 I = I . NR EG 

I F( NB( I ) .EO • 0 1 GU TO 60 

(J( I ) = I — G ( II ♦ AGA * All) )/ SO 

Till = 1 C I » .♦ CON * UN) 

IH TII). LE. 0.0 CALL CLEAR! I ) 
60 CONTINUE 
CALL DOSE 

CONST = ( 0 S— UiJ S ) / A A 
CL EC I = 1 .NREG 
I F ( NB( I J .60. 01 GO TO 80 
TII) = III) ♦ CONST * All) 

Iff TCIJ.LE. 0.0) CALL CLtARII) 
80 CONT INUE 


CALL WE IGF T 
CALL DOSE 

WRITE I6.6CC) IT. *1. DOS. (DID. 1=1. IMAX) 

write it.fcon tun. i=i,nregs 

I F ( ABSUT-b)/ U T -EPS) lie. 110, 120 
120 CONTINUE 

lie go io mi 

600 FORMAT! 1HC/7H0 II = I3.7X.5Hbf = I PE 1 2 . 4 , 7X ,6HCGS * 1PE12.4/ 
1 7HCi)(l) = 1P9E12.4/IYX.1P5E12.4)> 

601 FORMAT ( 7 EOT I I ) = IPSE 12.4/1 7X, IPSE 12.4) ) 

END 


1 I B FTC INPUT T DECK 

SIBROUIINE INPUT 


COMMON 

N K f: G « MAX. 

IMAX. DES. 

CON, CA. 

EPS. EPSD , f 

1 COS. AA. 

w r . gg * 

AG. Jf. 

NB 1 25 ) , 

N I J ( 25 ) « 

2 £125). 

AI 25). 

C ( 25) , 

E MU 1 25 .25) 

* 

3 T I 25 1 • 

RHOI 25) . 

G 1 2 5) * 

U 1 25 ) 


R E ACC 5.*GC) 

NREG. MAX. 

IMAX. DES. 

EPS. EPSD, 

CON, CA 

WRITE! 6. 500 ) NREG , MAX . 
READ I 5.4C1 ) (TII). I 

READ I 5.40 ) I RHOi 1 ) , 1 
Rf AC I 5.4C3 ) (N6I1 ). 1 

I MAX, DES. 
= 1 .NKEG ) 

=1 .NREG) 

=1 .NREG) 

EPS. EPSD, 

CCN , CA 


WR ITEI 6* 5 10 I ( l, [(I). KHC(l).NBU). 1 = 1. NR EG ) 

WR I TEI 6 . 550 ) 

CC 7 J = 1 *NKEG 

R F AC ( 5. 0 1 ) (EMU(l.J). 1 = 1.1 MAX ) 

7 WRITE 16.510 J.(EMUU.J). 1 = 1 .IMAX) 

READ I5.AU1) (U(I), 1=1. IMAX) 

K FAC (5. 40/ ) INIJ(i). 1=1. IMAX) 

C N1JII) IS THE KEG1CN NUMBER ASSOCIATED WITH THE ITH DOSE TERM 
Cf, S 1 = 1 , IMAX 
HE = 0.0 
DC 8 .1 = 1 .NKEG 
fl BF = KB * EMU I I, J)*I( J) 

S CII) = DU )*E XP IBB ) 

WR I TE 1 fc » 5 20 ) (I, C(l).. Dll). NIJUJ*. 1 = 1. IMAX) 

400 FCRMATI 3 15. 5F10.4) 

401 FCRMATI 7E1C.4I 

402 FCRMATI 14151 

403 FCRMATI 25 I ) ) 

500 FCRMAUSFl NKEG = 13/ SH MAX = 13/ 9H IMAX = 13/ 9H OES = 

1 1PF12.4/ SH EPS = 1PE12.4/ SH EPSD = 1PEI2.4/ 9H CON = 

2 1PE12.4/ SH CA = 1PE12.4) 

510 FCRMATI 1 F0/34HCKEGIUN Til) KHUII) NB I 1 ) / I I 7 ,2 F10. 3 , 17) ) 

530 FCRMATI IHC/24HC 1 CII) DID M J/ II 5 . 1 P2E 12 .4, I 5 ) ) 

550 FLRMAM ; 2ChCREGICN-J MUII.J) ) 

570 FCRMATI 16. JX. IP 9E 1 2 .4 / 1 SX , IP 9E 12 .4 ) ) 

R F T URN 
ENC 



TABLE IV. - Continued. OPEX-II LISTING AND OUTPUT FOR SAMPLE PROBLEM 


llBFIC INI1XX LECK 


L 

i 

3 


SLHPQbl INE 
Cl MM ON 
COS. A A . 

C( 25), 

It 25) . 
wR I TEC 6. 1 70) OES 
Q = 1 .0 

CC S5 K * 1 . MAX 


INI T 

NR EG * MAX. 
XT. CiG • 
AC 25) , 

RHOC 25) . 

DCS. 


IMAX, OES, 
AG. II. 

C C25) . 
GC25) • 

LT • C0C1). 


CON, CA, 
NB125) , 
EMU C25 .25 ) 
UI25I. 

1 = 1 . i M. A X ) 


EPS, EPSO, 
KI JC25), 


CC SA = COS 
A/A = SORTCAA) 

CCNST = CA*U /AAA 
IF ICOS-CtS) 62, 1C. fcC 
60 CCNST = -CONST 
62 CO 70 1=1. KUEG 

IM NHIIJ .Eg. 0) GC 70 7C 
TCI) = HI) ♦ CONST' * All) 

HI TCD.LE. 0.0) CALL CLEAR11) 

7C CONTINUE 
C/LL OUSE 

If( Art SC COS— UE S ) /OE S .LE. EPSO) GC TO 10 
If< C (f)0S-C2£S)«(U0SA-CES< ) * LE - 0.0 ) Q = 0.5 * 0 
S 5 CENT IN UE 
K = 777 
10 C/l.L WEIGH 
CALL OOSE 
WPITEC6, ICC) K 

100 FORMAT C iiHUlNIUAL IH1CKNESSES ELUNU AFTER I3.11H ITERATIONS) 
WPITEC6.U0) mi). 1=1, NREG ) 

1 6C FORMAT 1 2 1 FCCALCOLA IEU 1NI1IAL IHICKNESSES/ (1PSE12.4)) 
WR1TE16.17C) OES. UOS. WT , CO Cl) . 1=1, IMAX) 

1 70 FL'P MAT! 7FCCES = IPE12.4, 6X, 6hOCS = 1PE12.4, 6X. SHtaT = 1PE12.4/ 
1 7F0DC I >= 1PSE12.4/ C 7X, IPSE12.4) ) 

RFTURN 

END 


ilBFTC ILBANX DECK 


C 


TMAX. OES, CON. CA, 

AG, II, Nb C25) , 

C C 2 5) . E MU C 25 ,25 ) 

G ( 2 5 ) . UC25) 


SIERGUTINE CLEARCJ) 

COMMON NREG, MAX 

1 COS, AA. R T , GG. 

2 r. C 2 5 ) , A ( 2 5 ) . 

3 TC25J, KHCC25). 

TFt JTH REGION HAS JL'ST BEEN WIPED OUT 
TUI = 0 .C 

N Ft J ) = C 
DC 5 1=1, IMAX 

I F ( J «EU. MJt I ) ) CII)= C-0 
5 CC NT IMUE 


R f TURN 
ENC 


EPS, EPSO, 
M J C 25 ) , 



TABLE IV. - Continued. OPEX-H LISTING AND OUTPUT FOR SAMPLE PROBLEM 


iiejftc tdsexx ceck 


SlBKOUT INE 
CCMMUN 

1 rns# a a » 

2 C<25) • 

3 14 25). 

CCS - O.C 


LOSE 

NREG. FAX# 
WT. GG • 
At 25 )# 

RHO I 25 ) • 


IMAX, DES, 
AG# II# 

C ( 2 1) * 

G 1 2 5 ) # 


or 40 1 = J • I MA X 
HE = O.C 
CL 30 J= J • NK6G 
30 Ef * RB < EMU4 I # J>* It J) 
cm - = Cl I) * EXPI-BB) 

AO DCS = COS + 01 I ) 

DCS IS TOTAL DOSE 
AIK I * 4C-0GSE ) / (0— XK ) 

A t = A. A 


AC- = A.G 


A/ = 0.0 
AG - O.C 

DC- t 0 K = l.NkfcG 
AIK) = 0.0 

JH NKIKl .ED. 0) GO TO 60 
CO SO I = I • IMA X 

50 AIK) = AIK) - EMU ( l <K) * DU) 
A A = AA ♦ A t K )**2 
AG * AG ) A(K)*G(K) 

60 CCNTINUE 
R E TURM 
ENC 


CON# CA# EPS. EPSD. 
NB125). M J(25) i 

EMU125.25) , 

U( 25 ) 


i I RF 1 C USPHR CECK 


SlBROUTINE WEIGHT 

C CCG1NC, FOR OP EX SHIELDING OPTIMIZATION CCOE 

C SPHERICAL SHELL GEOMETRY -G.P. LAHTI -NASA-LEWIS 

f. FCK PLANE SLABS BUUNUEU BY CONE OF HALF ANGLE THETA, 

C THIS SUBROUTINE MAY BE USED SIMPLY BY REPLACING 

( THE TERM FOKP I BY 4*( IAN( THETA) )**2 

C/TA Ft) HP I / 12.56636/ 


C IMENS ICN 

ROUT (25 ) 



CCMMON 

NR EG « MAX, 

IMAX, UES# 

CON, CA. EPS, EPSD, 

1 COS, AA, 

w T • GG # 

AG# I I. 

Nb (25 ) , NIJI25), 

2 C( 25 ) . 

A I 2 5 ) • 

C (25) . 

EMU 125 ,25 ) , 

3 1(25). 

KHOl 25) . 

GI25). 

U 1 25 ) 


C CALCULATE GLTER RAD I LS OF EACH REGION. ANO WEIGHT 

K(tll)h Til) 

RPR= Rill T ( I )**3 
WE IGHI = RHO ( 1 )*RRR 
CC 5 . I =2 #NK EG 
R2 = RRR 

RCUTI I ) = HOI 1 I 1-1 ) ♦ I ( I ) 

K RK = ROL I I I )** 3 

5 WE IGF I ■= Wfc I GHT 4RHU ( 1 )* (KRK-H3) 

WI = WEI (FT*FORP 1/3. 

C CALCULATE PARTIAL WI DEK I WA TI VES 4DW/DXI) - GUI 
C GC = G.G 


CC S 1=1# NK EG 
RF = C.O 
GII) = O.C 

I F I NB< I) .TO. U ) GO TG S 
CC F J = I * NR EG 
H 2- RR 

Rh= KOLT ( J 2 

H GII) = GII) ♦ RHO I J)*IKK-R2) 
GII) = Oil) * FCRPI 
GC = GG -» 04 I )**2 
S CCNT IlML'E 
RETURN 
ENC 



TABLE IV. - Continued. OPEX-II LISTING AND OUTPUT FOR SAMPLE PROBLEM 


NR EG 

= 

11 

HAX 

= 

50 

IMAX 

3 

12 

CES 

= 

2.0000E 00 

EPS 

= 

l.QOCOE— 03 

EPS C 

* 

1 .00 00 E— 03 

CON 


7.CCC0E-01 

CA 


7. 00 ( OE— 0 1 


REGION Til) RhGt I ) 

1 26.000 9.951 

2 2.500 6.64 7 

3 0.600 16.763 

4 11. COO S.234 

5 17.5CC f . 750 

6 ?.U00 1S.3CC 

7 14.000 C.75C 

8 5.0GC 1S.JCG 

9 10.00C C.75C 

10 3.5CC IS. 300 

11 5S.5CC (.750 

P EG ICN- J MU 1 1 » J ) 


1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

c. 







2 

c. 

C. 

c. 

0. 

0. 

0. 

0. 

0. 

0. 


c. 

0. 

c. 







3 

c. 

c. 

c. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 

C. 

G. 







4 

c. 

c. 

C. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

c. 







5 

1.34 iOE-C 1 

1 .890CE— C2 

1.4 7C0E-02 

1.4700E-02 

2. 1200 E— 02 

2.4000E-02 

3.6290E— 01 

2.2970E-01 

2.01 80E-01 


1.8E50E-C1 

1.8370E-01 

1. 739CE-01 







6 

2.48CCE— Cl 

7.74COE-01 

7. 7200E-01 

7.7200E-01 

8.0500E-01 

7 . 9500E-0 1 

2.7190E-CI 

1 . 8120E-01 

2.3860E-01 


2.5 3606— 0 1 

2.5430E-01 

2.5770E-01 







7 

1.36SCE-01 

2.2flCCE-C2 

2.0400E-02 

2.04006-02 

1 .67006-02 

2.0600E-02 

2.1000E-02 

2.5000E-02 

2. 25 70E-01 


1.9830E-01 

1.9 ICO E— 01 

1. 795OE-01 







fl 

2.3 1 6 OE-Xl 

7 . 70CCE-01 

7. 66CGE— 0 1 

7.6600E-01 

7.9000 E-01 

7.8 lOOE-O l 

7.9800E— 01 

8.2700E-01 

1. 7790E— 01 


1 • 255CE-0 1 

2.0520E-C1 

2.3660E-01 







9 

1 . 2 180 6— C l 

2.4CCOE-02 

2 • 20GGE— 02 

2.2000E-02 

2.0100 E-02 

2.3900E-02 

1.8700E-02 

2.4500E-02 

2. 3000E-02 


2.5CCCE-C2 

2.231 OE-Oi 

i. doaoE-oi 







10 

2. 3060E-0 1 

7 • 6600E- Cl 

7. 6300E-0 1 

7. 63006-01 

8. 0500E-0 1 

7. 7600 E-01 

7.8200E-01 

8.1200E-01 

8. iOOOE-Ol 


6.51 COE-0 1 

2.407UE-01 

1. IG70E-01 







1 1 

1.18 70E-C1 

2.44CCE-02 

2 . 2600k— 02 

2. 26006-02 

1 .9600 E-02 

2. 5200E-02 

2.3200E-02 

2 • 0600 E— 02 

2.3200E-02 


2.9LCC6-02 

< . 3 IGCE— 02 

2.S7GCE-U2 








1 CC11 D1JI NI J 

L 2.942SE 05 2.43C0E-C2 1 

2 4.9449E C 3 2.03CUE-C3 1 

3 2.3SS7E Cl 1.S600E-C3 2 

4 2.6S36E (3 2.2CCOE-03 3 

5 3.5SFIE 05 2.04C0E-C1 4 

6 1.1522E C 4 5 . 04 COE— 03 4 

7 2.1S73E (h 9.21CUE-C2 6 

fl 2.25216 04 S.74C0E-C3 6 

9 S.5C79E 04 S.88C0E-02 £ 

10 2.C543E C4 2.78CUE-02 8 

11 1.1C3SE 05 2.01 COE- C 1 1C 

12 2.8426E C4 S.47COE-C2 1C 

CFS = 2.00CCE OC DOS = 7.6*676-01 U 1 = 3.5939E 07 

cm* 2.4 ioOt— 02 3.0300E-C3 l.SoOOE-03 2.2000E-03 2.0400E-01 5.0400E-03 9.2100E-02 9.7400E-03 9.8800E-02 

2.7H006-C2 2.01U06-C1 S.47CCE-02 

INITIAL 1H1CKNESSES FOUND AFTEk 11 UEkAlIGNS 

(ALC.IU.AIEC 1NI1 l AL THICKNESSES 

2.60CCE Cl 2.500CE CO 6 • CCCCE-C 1 1.1000E 01 1.77CIE 01 6.5022E 00 1.3850E 01 4.4536E 00 S.8800E 00 

2.H226E CC 5 .9466E 01 

CtS = 2.COGCE CC DCS = 2.COC36 CO H = 3.2192E 07 

Clll* 3.8S60E-C2 1 .1 5396— C2 7.40456-03 8.3112E-03 0.17O3E-O1 1.9&33E-02 2.9957E-01 3.0622E-02 2.2937E-01 

fc.455Sfc-C2 3.2S58E-CJ 1.4374E-C1 


N8( I 1 
0 
C 
C 
C 


17 



TABLE IV. - Continued. OPEX-O LISTING AND OUTPUT FOR SAMPLE PROBLEM 


IT - 

1 Vl 

= 3.15*26 

Cl 

COS 

= 2.00C5E 

00 





Cl I ) = 

4.19SIE-C2 
7 .48C 7E-C2 

9.6 44 St— 03 
3 .6601 E— C 1 

6.17C9E-03 
1. 545JE-01 

6.9265E-03 

6.7567 E— 0 1 

1 .63 1 IE— 02 

3.4357E- 01 

3.50 56E-02 

2.69386-01 

T 4 I ) — 

2 • 60C-GE C-l 
2 . 65S 56 CO 

2..50CUE CO 
5.S234E 01 

6 . COCUb-O 1 

1.1000E Cl 

1.7360E 01 

6.898 2E 00 

1.3516E 01 

4 .44226 00 

9.6920E 00 

JT * 

?. VT 

= 3.1772E 

C 7 

DOS 

= 2.0CC4E 

00 





Cl I » = 

4. 5. 1 63 E— 02 
£ . 1425fc-02 

7.57526-03 

4.0285E-C1 

5.C914E-03 
1.677 IE— 01 

5. 7148E-03 

5.5197E-01 

1 .34146-02 

3.8134E-01 

3 .09376- 02 

2.9864E-01 

1111 = 

2.600Gfc Cl 
2.66G2E CC 

2.5UC06 CO 
5.69206 01 

6.GUCCE-01 

1.1000E Cl 

1.7056E 01 

7.2748E 00 

1.3174E 01 

4.38296 00 

9.5066E 00 

IT = 

3 V T 

* 3.1652E 

C 7 

COS 

= 2.0004E 

00 





Cl I ) = 

4 .S506E-02 
H.C7 18E-02- 

6.654IE-C3 
4 .27856—01 

4.3631E-03 
1.83 35E-01 

4. 8974E-03 

4.67 72E— 01 

1. 1471E— 02 

4. 08 20 E— 01 

4.19016-02 

3.0352E-01 

T I I ) = 

2.60C0E Cl 
2.741:6 E CC 

2.5000E CO 
5 • 846 IE 01 

6 . CO COt - 01 

1.1UOOE Cl 

1 .6848 6 01 

7.5684E 00 

1.2842E 01 

4.23C9E 00 

9.34406 00 

IT = 

4 VT 

= 3.15536 

C7 

r 

ous 

=» 2.0003E 

CO 





Cl I ) = 

5.4284E-C2 
7 .4G48E— 02 

6 . 42 9 9fc — 03 
4.6407E-01 

4 • 0862E— 03 
1.55586-01 

4. 5665E-03 

4.3372E-01 

l • 0738E— 02 

4.19 39E— 01 

4.36006-02 

2.85806-01 

TU 1 = 

2.6000E Cl 
2. 5 34 46 00 

2.5000E CO 
5.7935E Cl 

6 .000 Cb— 01 

1.1000E 01 

1.68 0 3 E 01 

7.7148E 00 

1 .25666 01 

4 .0 1C 96 00 

9.24776 00 

IT = 

5 VT 

= 3.1463E 

C7 

cos 

= 2. 0003 E 

00 





01 I ) = 

5.9147 E— 02 
6. 8444 t— 02 

6.3 109E-03 
4.60C4E-CI 

4.C056E-C3 

2.12C5E-CI 

4. 4961 E— 03 

4.2181 Ef 01 

1.0526E-02 

4. 19 55 E— 01 

4.455 56-02 

2.69396-01 

Till* 

2.6000E Cl 
3.09?lf 00 

2.50CCE CO 
5.73346 Cl 

6 • 00 00b— 01 

1.1000E Cl 

1.6870E 01 

7 . 795 IE 00 

1.2329E 01 

3. 82 1 96 CO 

9.2176E 00 

IT = 

6 VI 

= 3.137SE 

C7 

DOS 

= 2.0003E 

00 





Cl I ) = 

6. 185 56 — 02 
6.515 78-02 

6. 2866E- C3 
4.EoSdE-01 

3.9862E -03 
2 . IV 5 IE- Cl 

4.474 3E—C3 

4 • 17 06E-0 1 

1.04 72E-02 

4. 1628E-0 1 

4.54886-02 

2.6057E-01 

T ( 1 ) = 

2.60CCE Cl 
3.188C6 00 

2 • 500 Cb CC 
5.67C4E 01 

6 .OOCOb-Ol 

1.1000E 01 

1.6994E 01 

7.8688E 00 

1.2120E 01 

3.6602E 00 

9.23906 00 

IT = 

7 V I 

= 3.129SE 

07 

oos 

= 2.0003 E 

00 





G( I ) = 

6.8410E-C2 
6.18 1fct-C2 

6 .3 166E-C3 
4.E 7C56-C1 

4 . 0029b— 03 
2.22S6E-01 

4 • 49 3 IE — 03 

4. 1657E-01 

1 .05 106—02 

4.12 73E-01 

4.6603E-02 

2.56756-01 

T 1 I ) = 

2.60CCE 01 
3.21726 00 

2.50CCE CO 
5.6063E Cl 

6 . 00 00b — 01 

I.IOCOE 01 

1.7147E 01 

7.9435E 00 

1.1940E 01 

3.56566 00 

9*29556 00 

IT = 

8 VT 

= 3 . 122 IE 

C7 

COS 

= 2.00C3E 

00 





Cl I 1 = 

7.2956E-C2 

6.33178-02 

6 .37646— C3 
4.E337E-C1 

4.03 63h— 03 
2 -23566—01 

4.5306E-03 

4.1809E-01 

1.0591E-02 

4.G941E-01 

4. 78466-02 

2.55816-01 

Tl I 1 = 

2. 60 COE Cl 
3.2583E OC 

2.5GCCE CC 
5.5416E Cl 

b.CUOOE-Ol 

1.1000E 01 

1.7311E 01 

8.0188E 00 

1.1784E 01 

3.4758E CO 

9.3693E 00 

II = 

5 VT 

= 3.1146E 

C7 

COS 

= 2.00C3E 

00 





Cl I ) = 

7 .7624E-02 
6 .3237E-C2 

6.446 IE— 03 
4.77736-01 

4.C77JE-03 

2.2366E-C1 

4.5765E-03 

4.2054E-0 1 

1.0692E-02 

4.0627E-01 

4.91476-02 

2.56316-01 

T 1 I ) = 

2 .60CCE Cl 
3.264 IE CC 

2.50CCE CC 
5.4 7 t 46 Cl 

6 • COOOE-Ol 

1. 10 00E 01 

1.7478E 01 

8.09 30 E 00 

1.1646E 01 

3.40356 00 

9.44976 00 


18 





TABLE IV 

. - Continued. 

OPEX-n LISTING AND OUTPUT FOR SAMPLE PROBLEM 


IT = 

10 V T 

= 3.1C73E 

C 7 

DCS 

= 2.0003 E 

00 





Cl I ) = 

8.2493E-C2 
6 .34C2E-C2 

6 . 52Q1E- 03 
4.711 3E— Oi 

4.12C9E-03 

2.2261E-0L 

4.6255E-03 

4.2335E-01 

1 .06 01E-02 

4.0332E— 01 

5.0475E-02 

2.5747E-01 

T ( I > = 

2.60CCE Cl 
3.2614E CC 

2.5QCCE CC 
5.41C7E 01 

b . OOOOt-C 1 

1.1000E 01 

1 . 764 IE 01 

8 • lfc 53 E 00 

1.1525E 01 

3.3365E OC 

9.5312E 00 

IT = 

11 t 7 

- 3.1C02E 

c? 

OCS 

= 2.CCC3E 

CO 





Cl t ) = 

C .761 2 E— C 2 
6 . 36 S <E— C2 

6. 5 94 8E- 03 
4.64C2E-U 

4.1650E-C3 
2.21 11E-01 

4.6750E-03 

4.2625E-01 

1 .09 1 IE-02 

4.0058E-01 

5.1818E- 02 

2.5889E-01 

11 1 1- 

2. 60 00 E Cl 
3.2538E 0C 

2-5UC0E 00 
5.3444E Cl 

6 . OOCOt-O 1 

1.1000E 01 

1 . 780 IE 01 

8.2355E 00 

1.1419E 01 

3.2765E 00 

9.61 12 E 00 

IT = 

12 v> T 

= 3.C933E 

C7 

OCS 

= 2.00C3E 

00 





Cl 11 = 

9.30C8E-02 
fc .405CE-C2 

6.6otioE-C3 
4 • 5.6CCE— C 1 

4.2084E-03 

2.1929E-01 

4. 7236E— 03 

4.291 IE-01 

1 • 1020E- 02 

3.9809E-01 

5.3174E-02 

2.6038E— 01 

71 I )- 

2 . 60C Ofc 01 
3.242 ?t CC 

2.500CE UC 
5.2778E Cl 

6. CUCCt-Cl 

l.iooot oi 

1.7954E 01 

6.3039E 00 

1.1327E 01 

3.2220E 00 

9.6886E 00 

IT = 

13 1* I 

= J.CEfc/E 

07 

DOS 

= 2.0003E 

00 





Cl I » = 

9 .€6S 7E-C2 
6 .445 IE-0'2 

6.7407E-C3 

4.46S5E-C1 

4.25C7E-U3 

2.1721E-01 

4.77 I2E-03 

4.3189E— 01 

1 .1126E-02 

3.95E5E-C1 

5 .4544E-02 

2.61 84E-01 

71 I ) = 

2.6000E Cl 
3.22tae cc 

2.50CCE CC 
5.21C7E Cl 

6. COCOt-ul 

1 . lOOut 01 

1 .8102E 01 

8.3704E 00 

1.12 50E 01 

2.1726E 00 

9.7628E 00 

IT - 

14 * 1 

= 3.CEC2E 

C 7 

COS 

= 2.0CC3E 

CO 





ci n = 

1 .04696-01 
6 .48e4t-C2 

6.eiC8E-Ci 
4 . 4 1 il E- C 1 

4.291 7E-03 
2 .146*6-01 

4 . 8 l 72fc -03 

4.3456E-01 

1 .12 30E-0 2 

3.9386E-01 

5.5931E-02 

2.6324E-01 

im-» 

2.60CCt Cl 
3.2124E CC 

2 . SUCGfc CC 
5.1433E 01 

6 . CUCOt-Ol 

L.loOOE 01 

1.8244E Oi 

8.4355E 00 

1.1187E 01 

3.1277E 00 

9.8334E 00 

n = 

15 fc r 

= 3.C74CE 

C 7 

ocs 

= 2.0QC3E 

00 





Cl I 1 = 

1 . 1099E-01 
6.534JE— 02 

6. 1 7fc fit- C3 
4 .33C6E— Cl 

4.3312E-03 

2.1234E-01 

4. 8616E-03 

4.3712E-U1 

1. 1331E-02 

3.9214E-01 

5.7 34QE- 02 

2.64576—01 

TUI* 

2.6UCCE Cl 
2.19351: OC 

2 • 50C0E CC 
5.C756E Cl 

6.C0C0E-01 

l . 1000E Oi 

1 • 838 IE 01 

8.4992E 00 

1.U37E 01 

3.0871E CO 

9. 9003 E .00 

IT = 

16 k 1 

= 3 . CfcCCE 

C7 

OCS 

= 2.00C3E 

CO 





U 1 ) = 

1 .17o lf-01 
6 .5823E-02 

6.S448E- C3 
4.24E 7fc— 0 1 

4.3o55E-u3 
2.C*5 4t-Cl 

4. 904oE-C3 

4.3957E-01 

1 .1429 E-02 

3. 9068E— 01 

5.8777E-02 

2.6580E-01 

7 < n = 

2 .ooooe ci 
3.1719E CC 

2 . 500 C£ CC 
5.C071E Cl 

6. CuCUE-Ui 

1 . 1O00E 01 

1 . 85 1 4t 01 

8.5620E 00 

1.1102E 01 

3.0504E CO 

9.9634E 00 

IT = 

17 H 

= 3.C622t 

Cl 

cos 

= 2 -0003E 

00 





c i n = 

1 .2** 5 3 1— 0 1 
6. 6325 E- 02 

7 .00 SUE- C3 
4 . 164 it— Cl 

4. 4U6oL-U3 
2 .C652E-C1 

4 . 94o 1 E— 03 

4.4193E-01 

1.1524E-02 

3.8947E-01 

6.02 50E-C2 

2.6694E-01 

T 1 I ) = 

2.6OC0E Cl 
3.147 it CC 

2.50CCE CC 
4.S396E Cl 

6 .coccc-oi 

1. 10C0E 01 

1.8o42E 0 1 

6.6240E 00 

1.1079E 01 

3.0176E 00 

1.0023E 01 
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TABLE IV. 

- Continued. 

OPEX-n LISTING AND OUT PUT FOR SAMPLE PROBLEM 1 


IT = 

18 UT 

= 3.05666 

07 

DOS 

= 2.0003E 

00 





C( I ) = 

1.31776-01 
6.68496— 02 

7.07146-03 
4. 07876-01 

4.44256-03 

2.0324E-C1 

4. 58656-03 

4 .442 IE— 01 

1 • 16 176— 02 

3.8852E-01 

6.1765E-02 

2.6798E— 01 

It I i= 

2 .60006 Cl 
3.12056 CC 

2.50006 00 
4-67136 .01 

6.CCC06-01 

1.10006 01 

1 . 87676 01 

8 .685 76 00 

1.10706 01 

2 .58826 00 

1- 0078E 01 

IT * 

19 k T 

= 3.05126 

07 

OOS 

= 2.0003E 

00 





C< H = 

1 .39316-01 
6. 7392 E— 02 

7 • 1323E— 03 
3.990 66- 01 

4.47 75E-03 
1 .55726—01 

5.02576-03 

4.4642E-01 

1.17G7E— 02 

3. 87836-01 

6 .33326—02 

2.6891E-01 

T ( I | = 

2.6OC0E Cl 
3.QSC46 (C 

2 -50CQE 00 
4.eo30E 01 

6.0CC0E— 01 

1.10006 01 

1.88896 01 

8.7473E 00 

1.1074E 01 

2.9622E 00 

1.0129E 01 

IT = 

20 k T 

= 3. 04 6 CE 

07 

DOS 

= 2.00C3E 

00 





0 { I 1 = 

1. 47146-01 
6.79‘4E-C2 

7.19196-C3 

3.90C4E-C1 

4.5116E-03 

1.5596E-01 

5 • 0640E-03 

4.4858E— 01 

1 .1795 E-02 

3. 87386-01 

6. 496 2E- 02 

2.6974E-01 

TI I 1 = 

2.60CCE 0.1 
2. €5716 0C 

2.50006 CO 
4.7346E 01 

6. COCQE-Ol 

1.1000E 01 

1.90086 01 

8.80926 00 

1.10926 01 

2 -9 39 IE 00 

1.0175E 01 

IT = 

; 1 k T 

= 3.04C56 

C 7 

DCS 

= 2.00C3E 

00 





Cl 1 1 = 

1.55266- (1 
6 .65346-02 

7. 25056-03 
3.8C61E-C1 

4.5451E-C3 
1.91 54E-01 

5.10 16E-03 

4.50716-01 

1. 18826-02 

3. 87 18E— 01 

6.6664E— 02 

2. 7045E-01 

II I » = 

2.60CGE Cl 
3 .02051 CC 

2 .500 06 00 
4.6662E 01 

6.CUC0E-01 

1.10006 01 

1.91266 01 

8.87176 00 

1 • 1122E 01 

2.9189E 00 

1 .02176 01 

IT = 

2? kT 

= 3.036 IE 

C7 

OOS 

a 2.00036 

00 , 





Cl 1 ) = 

1.63 636-01 
6. 91316-02 

7.3062E-C3 

3.71366-01 

4.57816-03 

1.87666-01 

5 s . 1386E-03 

4.52826-01 

1.19686-02 

3-8723 E-01 

6.8452E-02 

2.71 05 E-01 

im= 

2.60CCE (1 
2. 98026 CC 

2 • 500 OE 00 
4.55756 01 

6. 00006-01 

l. 10006 01 

1.9243E 01 

8 .93526 00 

1. 11656 01 

2.9012E 00 

L.0255E 01 

IT = 

23 VI 

= 3.C315E 

0 7 

COS 

= 2.0003E 

00 





C( I 1 = 

1 . 72 236- Oi 
6.5T44T-C2 

7.365 66-03 
3 .61686- C 1 

4.6L0UE-03 
1.631 46- Cl 

5. 17536-03 

4.54936-01 

1.20526-02 

3.87526-01 

7. 03406-02 

2.7152E-01 

1111 = 

2.60C06 01 
2 .53626 i)C 

2 .50006 CO 
4.52566 Cl 

6 . CO 006-01 

I.IOOOE 01 

1.9360E 01 

9.0000E 00 

1.12226 01 

2 .86596 00 

1.02876 01 

IT = 

24 UT 

= 3.02716 

07 

cos 

= 2.0003E 

00 





Cl I ) = 

1 .8104E-C 1 
7. 03 72 E- 02 

7.42276-03 
3.5179E— Cl 

4.6434E-03 

1.78366-01 

5.2120E-03 

4 . 5707 E— 0 1 

1.2137E-0 2 

3.88066-01 

7.2343E-02 

2 • 71 87E-01 

T I I > = 

2 .60 COE Cl 
2 .688 2E OC 

2.50CC6 00 
4.461 5E Cl 

fc.OOCOE-Ol 

I.IOOOE 01 

1.54776 01 

9.06666 00 

1.12916 01 

2.8726E 00 

1.03146 01 

IT = 

25 k T 

= 3.C2296 

C 7 

cos 

= 2.0003E 

00 





Cl I ) = 

1 .90 02 E— C 1 
7. 1C 1 76-02 

7.48006-03 
3 .4 1666— Ul 

4.67 d2E— 03 
1.73336-01 

5.24886-03 

4. 59266-01 

1. 22216-02 

3.88836-01 

7.4478E-02 

2.7209E-01 

Til 1 = 

2 .6000E Cl 
• 2.636CE CC 

2 • 500 CE CC 
4.39366 Cl 

6 . CU 006— Ul 

1.10006 01 

1.95956 01 

9.1353E 00 

1.13736 01 

2.8612E 00 

.1-03366 01 
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TABLE IV. - Concluded. OPEX-II LISTING AND OUTPUT FOR SAMPLE PROBLEM 


IT = 


l» T 

= 3. 0 i 6*36 

C? 

DOS 

= 2.00C3E 

00 





Cl I ) = 

1 .9912E-C 1 
1 . 1 hi it- 02 

? .537*6-03 
3.3 1316— Cl 

A. 70S*t-U3 
1 .6U06E-01 

5.2860E-03 

A. 61536-01 

1. 2306E-02 

3.b585E-01 

? «6766£— 02 

2.72176-01 

T ( 1 1 = 

<•600 Ut Cl 
/•7/S It 00 

2.S0CC6 CO 
<-325St Cl 

6.00006—01 

1.10C0E 01 

1.97156 01 

9.206AE 00 

1.14696 01 

2.8515E 00 

1.03536 01 

IT * 

21 

t 1 

* 3.015CE 

C 7 

DOS 

= 2.0003E 

00 





n i ) = 

2 . CH296- C 1 
7*2 34 lt-< 2 

7.59666-03 
3 .20756— CL 

4. 74326-03 
1.62556-01 

5. 32406-03 

A.6389E-01 

1 .239 IE— 02 

3.91 10E- 01 

7 .9222E-02 

2.72116-01 

T I I ) = 

2 -6ULC6 01 

2 .7 1 ?5fc CC 

2.5CCCE CO 
A. 25666 Cl 

6.00006-01 

1. 1000E Cl 

1.98386 01 

9.2805E 00 

1.15776 01 

2.8633E 00 

1.03636 01 

IT = 

26 ^ T 

* 3.0 11 At 

Cl 

DOS 

= 2.0003 6 

00 





DC I ) = 

2. 17A6E-L 1 
1.303 lb-02 

7. 656 86- 03 
3.CSS6E-C1 

A . 7 7 bOE-03 
1 .56 b2E-Cl 

5. 36316-03 

A. 66 38 6-01 

1.24 796-02 

3.92596-01 

8.18756-02 

2.71906-01 

111) = 

2 .60 CC E Cl 
2 . 6 506 E 00 

2.5CCCE CC 
A.1916E Cl 

6 • CO CUE— 0 1 

1. 1000E 01 

1.996A6 01 

9.35776 00 

1.16996 01 

2.83626 00 

1.03686 01 

IT = 

29 


= 3.0C7SE 

C 7 

COS 

= 2.00C3E 

00 





C< I ) = 

2 . 2t> 5 7£— C 1 
7-3731 £—02 

7.71 86E— 03 
2 .S9( IE- Cl 

6.ai3afc-03 

i.soeas-ci 

5. A03AE— 03 

A. 69026-01 

1.25686-02 

3.9430 E— 01 

8.67676-02 

2.71566-01 

Tl I 1 = 

2.60006 01 
2 . 57 6 7 E 00 

2.500CE CC 
A.1251E Cl 

e.cutuE-ci 

1.. 1000E 01 

2.0095E 01 

9.43 86 E 00 

1.1835E 01 

2.83036 00 

1.03666 01 

IT = 

20 

v r 

= . 3.CCA6E 

C 7 

OGS 

= 2.00C3E 

00 





Cl I ) = 

2.3554E— ( 1 
7 .44 ACE- C 2 

7 . 7829E-03 
2.6 76AE-C1 

A. 851 5E-03 
1 • AA 756—01 

5.AA56E-C3 

A.7186E-01 

l .26 616-0 2 

3.96236- 01 

8.78666-02 

2.71056-01 

Tl I ) = 

2.60C0E 01 
2.5011 t 00 

2.500CE CC 
4.C59CE Cl 

6 . CO CCt-Ol 

1.1000E Cl 

2.02306 01 

9.5233E 00 

1.19846 01 

2.82536 CO 

1.03576 01 

IT = 

21 

wi 

= 3 .001 At 

C7 

DCS 

= 2.0006 

00 





Cl I ) = 

2 . AA 2 fc E- C 1 
7.5171E-C2 

7 • 6A8SE- 03 
2. 7652 E— Cl 

A. 89036-03 
1 .38 A At — 01 

5. A891E-03 

A .7486 6—0 1 

1.27566-02 

3.98396-01 

9.12576-02 

2.70626-01 

Mil* 

2. 6000 t Cl 
2.A173E OC 

2.500LE CC 
3.9936E Cl 

6 . CO COE— 0 1 

1 . 10 00 E 01 

2,03716 01 

9.61256 00 

1.21466 01 

2.82126 00 

1.03626 01 

IT = 

22 

ul 

= 2 *S5tj 5t 

C 7 

DOS 

= 2.00C3E 

00 





C< I ) = 

2.52706-01 
7 . 5HS 76-0 2 

7.51SSE-C3 
2 ml 5C3E— 01 

A.932A6-03 

1.31596-01 

5. 53636-03 

A. 78226-01 

1.28586-02 

4.00726-01 

6.69666— 02 

2.69576-01 

M I ) = 

2.6O00E OJ 
2.32766 tC 

2.50CCE CC 
3 . S2E6E Cl 

6.COO06-01 

1.10006 01 

2.05196 01 

9.7056E 00 

1.23236 01 

2.81756 00 

1.0321E 01 


<C1» 6NUC5. tQF. REC- 00000 Ft L« 
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CONCLUDING REMARKS 


An optimization procedure for minimizing radiation shield weight has been described. 
The procedure, built around the steepest -descent code, OPEX-n, depends strongly on 
the validity of the dose -thickness relation. It has been observed that when care has been 
taken to accurately fit the dose -thickness relation, predictions of minimum weight con- 
figurations are quite good as verified by a necessary detailed proof calculation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 8, 1969, 

124-09-11-01-22. 
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Figured - Results of OPEX-II iterations. 
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